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We investigate relations between loop and tree amplitudes in quantum field theory that 
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involve putting on-shell some loop propagators. This generalizes the so-called Feynman 
tree theorems satisfied at 1-loop. Exploiting retarded boundary conditions, we give a 
generalization to £-loop expressing the loops as integrals over the on-shell phase space 
of exactly £ particles. We argue that the corresponding integrand for £ > 2 does not 
involve the forward limit of any physical tree amplitude, except in planar gauge theories. 
In that case we explicitly construct the relevant physical amplitude. Beyond the planar 
limit, abandoning direct integral representations, we propose that loops continue to be 
determined implicitly by the forward limit of physical connected trees, and we formulate 
a precise conjecture along this line. Finally, we set up technology to compute forward 

^ amplitudes in supersymmetric theories, in which specific simplifications occur. 
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1 Introduction 

Relations between loop and tree amplitudes are taking an increasingly important place 
in modern quantum field theory. To find such relations is actually a fairly old program 
[1]. Modern interest in it stems from the striking simplicity of on-shell tree amplitudes 
in gauge theories compared to the corresponding Feynman diagrams, which one aims 
to exploit in loop computations. The unitarity-based methods, pioneered in the mid- 
nineties in work by Bern, Dixon, Dunbar and Kosower [2], may be said to successfully 
meet this goal. These methods have now been successfully applied to multi-parton loop 
amplitudes in supersymmetric theories, but, also, in QCD. A non-exhaustive sample 
of references is [3, 4, 5, 6, 7]. 

The unitarity-based methods, generally speaking, exploit the fact that the discon- 
tinuities of a loop amplitude, or equivalently its unitarity cuts in all possible channels, 
suffice to determine it. The accuracy to which these cuts must be known may ulti- 
mately depend upon the theory; for instance, in QCD it is necessary to know the cuts 
in a neighborhood of D = 4 — 2e dimensions [8], while in supersymmetric theories, 4- 
dimensional data suffice [2] . Independent of these considerations, what is striking about 
these methods is that to find what a loop integrand was, one has to integrate it across 
unitarity cuts. In this paper we investigate possible more direct, purely algebraic, ways 
to determine the loop integrands. 

The idea is to try to determine loop integrands from their poles. Such poles are 
caused by loop propagators becoming on-shell, and so have the interpretation of on-shell 
particles circulating inside the loop. The corresponding data thus consists of on-shell 
tree scattering amplitudes in the forward limit, that is, zero momentum transfer for the 
particle in the loop. We will try to shed light on the physical origin of such relations, 
by trying to directly express loop amplitudes as phase space integrals over this data. 
Instead of unitarity, the relevant physical principle proves to be causality. 

The 1-loop version of this story was, in fact, worked out long ago by Feynman, and 
so it is sometimes known as "Feynman's tree theorem" [9]. It is useful to review it here. 
The first step is to use essentially contour integration in the frequency plane, to localize 
the frequency integrals to the poles of propagators. This succeeds in reducing any loop 
to on-shell phase-space integrals. The second step is to organize the combinatorics of 
the remaining part of the graphs, which are trees. Feynman could show at 1-loop that 
precisely the forward limit of physical tree amplitudes entered. 1 



1 Feynman used this result to exhibit the failure of longitudinal gauge modes to decouple, in gravity 
and Yang-Mills loops. To cancel them at 1-loop in an off-shell formulation he introduced particles 
now known as Faddcev-Popov ghosts. Somewhat amusingly, we are following the exact same road 
backward: we aim to move toward an on-shell formulation, where amplitudes are simpler, and, in 
particular, ghosts never appear. 



The contour integration step is actually subtle: it is not very naturally compatible 
with the familiar ie prescriptions. To see this, it is useful to follow Feynman and replace 
contour integration in the frequency plane by the one identity 2 

d 4 qG R (q)G R (q + Pl )---G R (q + p n ) q^ ■ ■ ■ q* = 0, (1.1) 



in which the G R are the retarded propagators: 
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In words, "a closed loop of retarded propagators is zero." This should be rather obvious 
in the time domain. At the same time, in frequency space, this follows from contour 
integration since the poles of all retarded propagators lie below the real q° axis allowing 
the q° contour to be deformed to infinity. This identity can be used to compute a loop 
with time-ordered propagators Gt, by substituting 

G R = Gt + G 

into it, where G K is an on-shell propagator accounting for the mismatch in the ie 
prescriptions (for precision, the "<" superscript refers to its support being on negative 
frequencies). Schematically, for a loop containing m propagators this leads to 

/ d D q {G T ) m = ! d D q [{G T + G<) m - {G T ) m ] . 

Expanding out the bracket gives a series of terms each one of which contains at least 
one on-shell propagator. After summing over Feynman diagrams, this can be organized 
in terms of physical data. 

We will not use this form of the tree theorem. This form contains (2 m — 1) terms 
on the right-hand side (not all of which may be kinematically allowed, of course), 
which can be a lot. The bulk of them comes from just changing G R s to Grr's and 
can be interpreted as just changing boundary conditions. Rather, we will avoid such 
contributions by considering observables with retarded boundary conditions. One can 
then hope that a simple analytic continuation will relate these observables to time- 
ordered ones, and this will indeed be the case. In any case, since the identity at play 
is due to causality, it makes perfect sense to try it on causal observables. 

What should we use as "causal observables"? Since retarded boundary conditions 
are a preferred choice in classical field theory, we turn to there for inspiration. In that 
context, a natural question is, under the influence of small perturbation, how do fields 
change in the future? Response functions answer this question, if we define an n point 
response function to be the solution to the field equations with retarded boundary 
conditions and (n — 1) infinitesimal sources, the n-th leg becoming the observed field. 
Perturbatively, classical response functions are computed just by summing over all tree 



2 We ignore any issue with convergence and arcs at infinity, which are absent in dimensional regu- 
larization. 
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diagrams as usual, but with the replacement Gt — > Gr for every propagator. The time 
flow in every retarded propagator, evidently, is such as to point toward the observed 
field. In configuration space a response function vanishes unless all (n— 1) sources are in 
the past lightcone of the n-th leg. Since the same question can be asked in a quantum 
theory, response functions can also be defined at the quantum level. Response functions 
may be familiar to some reader from many-body or finite-temperature contexts [10]; 
for the benefit of other readers, their definition and key properties will be reviewed 
below. For the moment, the most important aspect to keep in mind is their space-time 
support. 

For a 1-loop response function the tree theorem takes the particularly concise form: 

i£ loop (pi, • • • ,Pn) = EH 2 " f^ + ( q 2 +ml) rT + c 2 (pi, ...,Pn\ <A -q~ h ), (1-2) 

, 2. d D q 
q o " (2tt)» 

is the ordinary loop integral in D = 4— 2e dimensions and the sum is over the helicities 
h (or spin) of the particles. This form was obtained relatively recently by Catani et 
al. [11] using contour integration arguments. Morally, the integrand is to be evaluated 
in D dimensions, if one is to work within a consistently regularized theory. Both sides 
of the equation involve retarded amplitudes, for which there is one special "future- 
directed" leg (the same on both sides). We do not indicate this leg explicitly and the 
semi-colon refers to the fact that the last two legs on the right are amputated. For a 
loop diagram with m propagators, this form contains exactly m terms. 

It is actually easy to understand why this form contains only m terms. This is due to 
the boundary conditions. On-shell propagators have no notion of time flow so a product 
of two subgraphs separated by on-shell propagators could not possibly obey retarded 
boundary conditions. There must always be a connected tree of retarded propagators. 
This argument is very powerful, and it applies to any loop order: at £-loop, at most 
£ propagators can get cut in a causality-preserving fashion. On the other hand, using 
contour integration in £ frequency variables, one can always cut at least £ propagators. 
The conclusion is that, any £-loop amplitude can be represented as a phase space 
integral over exactly £ on-shell particles, times connected trees! Investigating such 
formulas is the point of this paper. 

At 1-loop (1.2) gives rise to a quite elegant picture of the quantum vacuum, as 
being filled with on-shell particles, whose interactions with us give rise to quantum 
corrections. We cannot see them, because they uniformly fill space-time, like a bosonic 
version of the "Dirac sea." But although invisible in that sense, they do generate tan- 
gible quantum effects. In this paper we will see quite explicitly that such a picture fails 
at 2-loop in general, a physical interpretation of the connected trees being unavailable. 
An exception will be for planar gauge theories, where introducing a very specific set of 
color-correlations for the vacuum particles will restore this interpretation. This nicely 
fits within the idea that planar theories are "classical" , albeit in an unusual sense. 

This paper is organized as follows. In section 2 we attempt to "guess" the higher- 
loop extension of (1.2) based on causality and unitarity. A clash between those two 



principles will be argued to obstruct physical integral representations. Specifically, we 
will find a unique integral representation but it will not involve a physical combination 
of tree diagrams. This explains why a 2- loop extension of (1.2) has never appeared 
in the literature. However, we will see how, in the planar limit, the unphysical sum 
over trees becomes equivalent to a physical one. In this paper we also make a proposal 
regarding non-planar theories. Abandoning Minkowski-space phase-space integrals in 
theses case, we propose that the forward limits of physical connected trees continue to 
implicitly determine loop integrals. For the reader's convenience, a precise statement 
in the form of a conjecture is given at the end of section 2. 

Section 3 contains a pedagogical introduction to response functions and causal 
boundary conditions in quantum field theory. Their relation to time-ordered functions 
is also explained. An efficient formalism to compute them, the Schwinger-Keldysh 
formalism, is introduced. Quite pleasingly, we will see that the Feynman rules of this 
formalism, when expressed in the appropriate field basis, directly generate the 1-loop 
formula (1.2) without requiring any work. 

In section 4 we apply this formalism to 2- and 3-loop diagrams. A certain special 
vertex entering the Schwinger-Keldysh formalism will enter, which will turn out to be 
related to the obstruction found in section 2. The evaluation of these diagrams will 
establish the 2-loop formula deduced there involving unphysical trees. Starting from 
this formula, we will prove the physical one valid in planar gauge theories. This section 
concludes with formal manipulations supporting our "uniqueness conjecture", which 
exploit complexified Minkowski space. 

In section 5 we work out a simple 2-loop example of the formalism: the 2-loop 
vacuum polarization in QED (or QCD). We calculate it starting from physical on-shell 
tree amplitudes. 

Finally, some technology to efficiently compute the forward amplitudes in super- 
symmetric theories is discussed in section 6. In this case the definition of the forward 
limit of amplitudes for S-matrix elements, as opposed to correlation functions, is fully 
worked out. This is simplified by supersymmetry cancellations. Another simplifica- 
tion is the fast decay of amplitudes at large loop momenta, which allows the use of 
factorization or recursion formulas. A few concrete examples are given in the simplest 
case of planar M = 4 supersymmetric Yang-Mills, which is clearly singled out as the 
simplest theory to which to apply this formalism. 

The tree theorems, as well as nontrivial implications of causality at tree-level, have 
appeared in the recent works [11, 12, 13]. 

2 Causality and unitarity 

We have just argued that causality implies the existence of on-shell integral represen- 
tations of loop amplitudes. We now wish to argue that by combining this information 
with unitarity, the form of these relations is uniquely determined. 

At 1-loop, (1.2) is pretty much the only answer one could get: one has to open 
exactly one propagator in the loop. For a 1-loop n-point amplitudes, this leaves out a 
(n + 2)-point tree amplitude as in (1.2). Assuming that combinatorics work out nicely, 
this tree amplitude has better be the forward limit of a physical tree amplitude. We 



will check below that this is the only possibility consistent with unitarity, but let us 
move on to 2-loop. 

A naive 2-loop exponentiation of the pattern would yield 
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^2-loop, classical^ pn) = 2- J[ ^T(-) 2 ^ f **+(??+<) 

x i2a(pi,...,p B ;rf 1 ,-gr fcl ,^ a ,-?2" ha )- (2-3) 

This turns out not to be correct, meaning that we have to think harder. 

The reason this is not correct is the following. This formula amounts to treating 
the vacuum as a collection of classical fluctuations: in it the fields obey the classical 
equations of motion, but they are averaged over a distribution of initial conditions. The 
distribution that exponentiates to (2.3) is a Gaussian, and there can be nothing wrong 
with that if the initial conditions are set at some early time where the interactions 
are adiabatically shut off. Furthermore, the two-point function read off from (1.2), 
(— ) 2h 7i5 + (q 2 ), gives each oscillator mode an average energy (—) 2h ^hu which would 
seem to be one's best guess at a classical description. The problem, now, is that (2.3) 
does not address any of the well-known problems of the classical theory. For instance, 
the energy distribution h^ is not dynamically stable - it would attempt to thermalize 
- making this expansion not even self-consistent. A second profound problem is energy 
positivity: in a classical theory endowed with "vacuum" fluctuations, energy would 
not remain in the "vacuum" but would transfer itself to any particle present on top 
of it. This, of course, does not happen for zero-point energy in quantum field theory, 
it being impossible to extract energy from a state of lowest energy, in agreement with 
observations. 

Since energy positivity is closely related to unitarity (it being the familiar statement 
that energies along a unitarity cut are positive), this second problem implies that (2.3) 
violates unitarity in a very bad way. Let us expound this in more detail. Consider 
the example of an electron's self-energy in QED. According to the classical model 
(2.3), the electron would acquire a large collisional imaginary part through inverse 
Compton scattering against classical fluctuations of the electromagnetic background. 
Diagrammatically this is the 2-loop graph shown in Fig. 1 (a). Note that this is 
2-loop graph, not contradicting the validity of a classical picture at 1-loop. Taking 
the imaginary part of this self-energy (which for retarded boundary conditions is also 
the lifetime of the particle) gives, for the term in (2.3) corresponding to (b) and (c), 
graphs in which the electron propagator is cut so this gives a unitarity cut. Now if 
both (b) and (c) are included, as dictated by (2.3), there will be a contribution in 
which the two photon energies have opposite signs in the loop, corresponding to an 
absorption-emission process, that is, inverse Compton scattering. This produces the 
linearly divergent decay rate 3 

^uv- 
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Energy from the classical "vacuum" is transferred to the electron at this rate. Evidently, 



3 The complete evaluation of (2.3) for Fig. 1 (a) also includes the vacuum fluctuations of the 
electron's field. This modifies the coefficient but does not remove the linear divergence. 
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Figure 1: (a) A 2-loop electron self-energy diagram. The photon energies flowing in 
it can be parallel or anti-parallel as in (b) and (c), each case having to be considered 
separately in order to respect energy positivity. In some other diagrams (d), the en- 
ergies in two cut propagators cannot be compared as they do not add nor subtract in 
any denominator. 



such a huge scattering rate would destroy any atomic system almost instantly; (2.3) 
fails at even the most basic qualitative level. 

Ordinarily, Feynman's ie prescription guarantees the absence of such unphysical 
contributions, which is why one never worries about them. What happens here is that, 
to cut-open the loop, we had to use causality which clashes with energy positivity (and 
with Feynman's ie prescription). The price to pay for working with causality-friendly 
propagators is that the unphysical region no longer cancels in individual diagrams. 
It must, of course, cancel in the end through some other mechanism. After all, the 
underlying quantum field theory is the same and we are merely asking a different sort 
of question. 

The only way we can correct this, considering Minkowski-space integral representa- 
tions, is to correlate the signs of the energies q\ and q 2 ■ This forces us to use a nonstan- 
dard weighting of tree diagrams so as to avoid the automatic symmetry (g, -B- — <&) for 
each of the <&. The most general weighting can be described as follows. For any tree 
diagram with four special external legs, (q±, —q±, q 2 , —q 2 ), there are two unique oriented 
paths which link q 1 to —q\ and q 2 to — q 2 (see Fig. 1). If they share a propagator with 
the same orientation (that is if some propagator carries momentum (qi+q 2 +p) for some 
p), we say that qi and q 2 are parallel, and if there is minus sign, we say that they are 
anti-parallel. It may also happen that they do not share any propagator, in which case 
we say that q\ and q 2 to be non-comparable. One linear combination of these cases 



must respect unitarity. Let us state it, and then explain it: 

- -=i L hi J * 
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connected trees 

1 + 1/3, qi and q 2 parallel, 
x ^ 1 — 1/3, qi and q 2 anti-parallel, (2.4) 

1, gi and q 2 non-comparable. 

The curious factors | ensure unitarity. To see this, consider a three-particle unitar- 
ity cut. Essentially, | is because only two energies in the unitarity cut can be controlled 
in a given term. Let us focus on the unphysical region in which two propagators have 
positive energy and the third has negative energy. In this region, for one term in (2.4), 
<7i and q 2 will have positive energy making it a "parallel" configuration. In two other 
terms, q\ and q 2 will have opposite energies making them "anti-parallel". Accounting 
for the sign flip for the imaginary part of a retarded propagator at positive and negative 
energies, the total imaginary part from this region will be 

— (parallel) + 2 (anti-parallel). (2.5) 

Unitarity (or energy positivity) requires that this be zero, which fixes the ratio of the 
first two terms in (2.4). The unitarity cuts in the all-energies-positive region (which 
is the same for a retarded and a time-ordered function, in the absence of a fourth 
propagator) determine the remaining coefficients; note that the average of the coeffi- 
cients just matches the classical expectation. The conclusion is that (2.4) is the only 
Minkowski-signature integral representation consistent with causality and unitarity at 
2-loop. 

Related problems with a classical model at 2-loop were also found by [14]. 

At 1-loop, the same computation gives a "1" instead of a "2" in (2.5), showing that 
(1.2) duly respects unitarity. 

Formula (2.4) will be proved in section 4 using a method not relying on unitarity. 
But this is not yet the formula we are looking for: its integrand is unphysical - we find 
that longitudinal polarizations and ghosts do not even cancel in it. 4 

2.1 Planar gauge theories 

The situation proves to be much better in planar gauge theories. This is for purely 
kinematics reasons: the color flow in a planar tree diagram is so constraining that, 
fixed the external momenta and color structure, the momentum flow in all channels is 
completely fixed. Hence one can replace the condition on the momentum flow in (2.4) 
by one on the external color structure. Equivalently, the problem with the instability 
of the classical vacuum is solved by a universal distribution of color correlations, one 
which can presumably be computed loop order per loop order. 



4 They do, of course, cancel when the integral is taken, for gauge-invariant observables, as usual. 




Figure 2: Color structures for forward amplitudes are in one-to-one correspondence 
with cut planar diagrams in open-string theory. 

Specifically, consider a gauge theory with a double-line notation, such as U(N) 
gauge theory, and consider the large N limit with all fields in the adjoint representation. 
Tree amplitudes can be (uniquely) expanded as 



-Rn(l> 
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J2 2# n - 2 Tr [£ ai £<M2) . . . £Mn)] R n {\ i a (2), ..., a(n)). 
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Extra legs added to compute planar loops are color-adjacent. At 1-loop, for in- 
stance, they can be inserted in n ways, R n+2 (. ■ ■ ,Qi,Qi, —Qi, Qi+i, • • •)• At 2-loop there 
are exactly n(n + l)/2 possibilities. They fall into two distinct classes: two disjoint 
pairs -R n+4 (. . . ,q 1 ,-q 1 , . . . , q 2 , -q 2 , ■ • •) or two nested pairs R n+A (. . . , q ly q 2 , -q 2 , -qi). 
It is easy to see that the first case includes all Feynman diagrams where q\ and q 2 are 
parallel, while the nested case includes all diagrams in which they are anti-parallel (the 
two cases get exchanged, evidently, when the sign of qi or q 2 get exchanged). Hence, 
following the above considerations, an Ansatz can be constructed consistent with both 
causality and unitarity: 
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4/3, sgn q\ = sgn g§, 
2/3, sgn ql = -sgn g°. 



(2.7) 



This formula turns out to be, formally, an exact identity. In section 4.3 we will deduce 
it from (2.4) and discuss some subtleties it involves. 

A useful mnemonic for (2.7) is that it is precisely (2.4) applied to planar open-string 
theory diagrams, as in Fig. 2; a similar version is obtained at 3-loop following section 
4.2. This is quite pleasing: it simply means that the (space-time) Schwinger-Keldysh 
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formalism, which will be used below to derive (2.4), is applicable directly at the level 
of individual open-string diagrams. 

2.2 Non-planar theories 

Equation (2.4), viewed as the sum of a classical and a quantum part, gives control 
over the quantum part. We can try to determine to what extent it is sensitive to 
unphysical data. In section (4.4) we will argue that its integral (the part with the |'s) 
is in fact fully determined by physical 3-particle cuts (only, at complexified momenta!) 
This is not unreasonable, since the only reason this term had to be introduced was 
to fix 3-particle cuts. This term vanishes for non-overlapping loops, which have no 
3-particle cuts. This leads us to make the following suggestion: the quantum part may 
be implicitly determined by the classical part. Lacking a good physical interpretation 
of such a result, we will be content with formulating a conjecture. 

This is best stated as a uniqueness conjecture. Specifically, suppose one has a 
representation 

<" l00P = E C n X J n~ l0 ° P ' (^ 

i 

as a sum over ordinary £-loop Feynman integrals I n (e.g. any product of 1/(5^ 9i + 
J^p) 2 -type denominators together with numerator factors), which are to be computed 
in Euclidean space as usual and then continued to give the retarded amplitude. 5 The 
coefficients C l n are rational functions of the external momenta only. This is the same 
sort of expansion as is familiar from unitarity-based methods, where one would take 
the J n 's to span a basis of scalar integrals such as tadpoles, bubbles, triangles, boxes at 
1-loop. For the present discussion it will not matter what the £-loop basis is; it could 
be complete or over-complete. Since both sides are Feynman integrals one can take 
their classical "^-particle cuts", defined as the integrand to (2.3), and compare them: 

E n 2hi+ - +2he Rn e + e M\ -QT h \ ■ ■ ■ , Qe\ -<u hi ) 

hi,...,he 

= E C l n x [classical ^-particle cuts] (/^~ loop (gi, —qi, . . . , q e , -qi)) ■ (2.9) 

i 

In general, (2.8) does not imply (2.9) since the right-hand side of (2.8) is defined only 
modulo terms which integrate to zero. Our uniqueness conjecture (supported in section 
4.4 for £ < 3) is the other implication: 

Conjecture. If (2.9) is an equality, then (2.8) is an equality. (2-10) 

The left-hand side of (2.9) is given by the forward limit of physical tree diagrams, and 
can thus, in principle, be computed directly. That means that if an Ansatz (2.8) for 
an integrand agrees with physical trees on all single-cuts, then its integral gives the 



5 Actually the same Euclidean amplitude, that is the same coefficients C l n , will determine both 
the retarded and time-ordered functions. The relevant analytic continuation is described in the next 
section. 
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correct loop integral. This gives a set of purely algebraic constraints to determine the 
basis coefficients C l n . 

Some clarifications are in order. First, any loop integrand is at best defined only 
modulo shifts in loop momenta (except in a planar theory, where origins can be defined, 
but since we are now considering non-planar theories this cannot be invoked), and the 
right-hand side of (2.8) is subject to this freedom. The cut propagators, however, 
provide a natural origin for the £ loop momenta and the comparison in (2.9) does not 
involve any choice of origin (symmetrization under permutations of the loop momenta 
is implied). 

Second, this allows the determination of loop integrands modulo integrands with 
vanishing ^-particle cuts. At 1-loop, we believe that any such integrand is either a 
polynomial or a massless bubble. 6 More generally, the classification of such integrands 
is not known but we did not need it. The conjecture is that any such integrand 
integrates to zero. 

Third, there exists integrands which integrate to zero but which have non- vanishing 
single-cuts. We will see an explicit example in section 6.3.2; in general, a class of 
examples are integrals which vanish upon symmetrical integration. Because there is no 
integration in (2.9), such terms do not necessarily vanish in it. However it is also clear 
that one could perform e.g. symmetrical integration at the level of (2.9), if one wished. 
Thus one is free to keep or ignore such terms as one wishes. To put it differently, 
what is to be judged "equal" in (2.9) depends on what one wants to achieve. Similarly, 
whether the basis in (2.8) has to be a full algebraic basis, as in the OPP method [15], 
or a more reduced basis, as in the conventional unitarity-based methods, depends on 
the amount of vanishing integrals one wishes to keep. 

Finally, it is important to note that matching all ^-particle cuts on the two sides, 
automatically implies that all unitarity cuts match. It is certainly always possible to 
match all ^-particle cuts. Thus the conjecture would be proved, if a proof that (D- 
dimensional) unitarity cuts suffice to determine any £-loop amplitude, existed for any 
theory. 

3 Response functions and Schwinger-Keldysh formalism 

We review the basic properties of retarded response functions and their computation 
within the Schwinger-Keldysh formalism. Our presentation follows a semi-classical 
expansion which may be more familiar from quantum mechanics. 

The response functions are important for the present considerations, as they dra- 
matically simplify the physics and help keep equations such as (2.4) intelligible. The 
reader familiar with them could skip directly to the next section; some useful references 
for this section are [10, 16, 17]. 



6 The cuts of massless bubble (a loop integral with two internal propagator with no mass and a 
massless external momenta) are equivalent to l/(q+p) 2 + 1/(<Z — p) 2 = ^-/2q-p + 1/ — 2q-p = because 
arbitrary shifts in the integration momenta are allowed. This is never a problem in practice, because 
massless bubbles integrate to zero in dimensional regularization. 
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3.1 Response functions 

Response functions describe the "response" to a perturbation. One considers a system 
initially in its vacuum state |0) and add a perturbation H — > Ti + f(x)0(x) to its 
Hamiltonian density. At time t the system will be in the state 

|^(t)) / = Te-i/^x/WOWio^ (3.H) 

where T is the time-ordering symbol and 0(x) is in the Heisenberg representation, as 
usual in field theory. By analogy with classical response functions, quantum response 
functions should describe the expectation value of an operator 0\{x) as a function of 
/, or, more precisely, as a series expansion in powers of /. This motivates the definition 
for the n-point response function R n (l, 2, . . . n): 



nrn-l 



R n (x 1 ,x 2 ,...x n ) = )...Sf( x y ^W" ^ 1 '^®'} 



(3.12) 
/=o 



It is possible to expand this definition in terms of vacuum expectation values of operator 
products. The 2-point function becomes the usual retarded product 

R2(x 1 ,x 2 ) = -i9{x\-xl) (0\[O 1 (x 1 ),O 2 (x 2 )]\0), 

while n-point functions are expressed in terms of series of nested commutators [10]. 
We shall not need such formulas. 

Important properties of these objects are: 

• They are Lorentz-covariant despite the explicit appearance of time in the 9- 
functions. For the two-point function this may be seen from the vanishing of 
commutators at space-like separated points. More generally, this follows from 
the definition (3.12): imposing the vacuum state at past infinity is a Lorentz- 
invariant boundary condition. 

• In momentum space the response functions are a simple analytic continuation of 
the Euclidean amplitudes. In this continuation, a positive imaginary part is added 
to all incoming momenta p 2 . . .p n . (The relative magnitudes of the imaginary 
parts is not important, because all discontinuities lie on the real axis.) This is 
discussed in details in [18]. In particular, in vacuum, response functions are an 
analytic continuation of time-ordered functions, since the latter are also obtained 
from Euclidean ones by a Wick rotation. (The relation between response and 
Euclidean functions does extend to general, finite temperature, contexts [18]. In 
that context retarded functions become the simplest objects, with time-ordered 
functions becoming neither particularly simple nor useful.) For e.g. a simple 
power law the relationship is just 

[p — te) I time-ordered """ ^P ~ te P > I retarded' 

In the general case, the analytic continuation represented in Fig. 3 should be 
used. An explicit example will be given in the next subsection. In general, it 
establishes that retarded and time-ordered amplitudes have the same information 
content. 
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Figure 3: The three-step analytic continuation which takes retarded amplitudes to 
time-ordered ones: from Minkowski (real q°) to Euclidean signature, then in Euclidean 
signature from the positive to the negative imaginary axis, then from Euclidean to 
Minkowski signature through a Wick rotation. 



3.2 Example: a scalar triangle 

To help reassure the reader that the integrals we discuss are indeed closely related to 
familiar loop amplitudes, let us immediately describe a simple but nontrivial example: 



(3.13) 



q 2 (q+Pi) 2 (q-p2) 2 ' 



Here we consider the case that p\ = p\ = for simplicity. The usual result, with 
time-ordered boundary conditions, is 

h 



-M-iz)- 1 -' 



(3.14) 



where cr is a function of e defined in Appendix B and p 3 = —p\—p2- Let us verify that 
(1.2) indeed produces a simple analytic continuation of that answer. 

According to the latter formula, the triangle integral should be equivalent to 



7T 



5(q 2 



1 



1 



1 



(2q-p 1 +p 2 1 ){-2q-p 2 +p 2 ) (2q-p 2 +pl)(-2q-p 3 +p 2 3 ) {2q-p 3 +pf){-2q-p l +p 2 ) 

(3-15) 
This sum is depicted in Fig. 4 (b). The ie prescriptions are also specified by the formula. 
With retarded boundary conditions one must choose one special leg, corresponding to 
the "advanced" one, which we will take to be p 3 . Then in all terms, time is to flow 
from the cut propagator toward p%. For instance, the first term is to be read as 



(2q-p ± +pl - ie(p?+g°))(-2g-p 2 +p| - Hp 2 -q )) ' 

(A convenient mnemonic device, which works for all terms, is to generate the ie's by 
adding small positive imaginary parts to p° an d p®-) As a rule, we omit the explicit 
writing of ie factors. 
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It is important to note that the ie's are only needed when a retarded propagator 
gets close to its mass shell. The sign of its imaginary part, opposite to the sign of its 
energy, is then a well-defined Lorentz-invariant quantity. 

The integrals are easily computed in a frame in which p\ and p 2 are within the z — t 
plane, so that 2p r p 2 = —p\p 2 with pf = p® ±pf. We begin with the second term in 
the bracket. Upon using the S- function to perform the q± integration in D dimensions 
it is reduced to 

dq + dq~8(q + q~) 1 



(47r) D / 2 r(l - e) J 2(q + q-y q+p 2 (q-p+ + q+p 2 + p\p 2 + ie(q° +p\+p%)) ' 

(3.16) 
Three possible cases must be distinguished. 

When both p± and p 2 are positive the pole occurs at negative q + ,q~- The sign 
of the ie prescription is dominated by that of {pi+p 2 )° > (in the rest frame of p 3 , 
q° carries half the energy), so the pole is below the real axis. Thus the integral for 
negative q + , q~ may be evaluated by rotating the two variables by 180° clockwise, while 
that for positive q + ,q~ is done directly, yielding for (3.16) 

r(l - e)r(2e) 2 1 + e 2 ™ 

e{Aix) D ^pl X |P3 ' 2 • 

When both p\ and p 2 are negative one finds the complex conjugate 

r(l-e)r(2e) , l + e-a* 
e(47r)^/ 2 pl ^ 2 • 

When p^ and p 2 have opposite sign, that is p\ is spacelike, the sign of the ie addition 
depends on whether it is q + or q~ which is larger at the pole, which is determined by 
the sign of q + ,q~- For instance if p^ > then for one sign the integral can be done 
by rotating q + clockwise and for the other sign it can be done by rotating q~ counter- 
clockwise. The result is therefore 



TV IX 



r(l-e)r(2e) e ™ + e- 

e(47r) D /y 3 lPsl 2 

The third term in the bracket (3.15) gives exactly the same result and the first term 
gives zero, there being no scale in it. To obtain the final answer the above results must 
thus simply be multiplied by 2. The three different cases nicely combine into a single 
formula, 

^(pg + iepg)- 1 - 

e z 

which is precisely the claimed simple analytic continuation of (3.14). 

In practice, computing any one of these cases is enough, as all three are related 
through well-defined analytic continuations. 

If we had used a different retarded boundary condition so that, for instance, p\ or 
P2 were the advanced leg instead of p$, the same computation would have landed us in 
a different branch of the result 



— x(p 2 3 - itpl) 
e z 
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again in complete agreement with the analysis of the preceding subsection. 

In summary, regardless of how results are obtained, it is trivial to convert between 
time-ordered and retarded boundary conditions. In the Euclidean region (spacelike 
invariants) there is in fact no continuation at all to be done. 

The reader interested in an example with a less trivial analytic form than (3.14) is 
referred to Appendix B. There we compute the triangle for general external momenta, 
which in D dimensions is expressed in terms of hypergeometric functions. Full agree- 
ment is again obtained, which follows from a non-trivial identity between Feynman 
parameter integrals of different dimensionality. 

3.3 Semi-classical limit using the Schwinger-Keldysh formalism 

Response functions admit a path-integral representation known as the Schwinger- 
Keldysh path-integral. This formalism is employed to study finite-temperature sys- 
tems, and even out-of-equilibrium systems. Here, we use it to study a theory in its 
vacuum state. But since the formalism is the same, it is worth noting that much in 
this section (and paper) should generalize to these other cases. 

This path integral runs over a doubled set of fields. It would be impossible not to 
double the fields to compute (3.12): one set is needed to evolve the ket, and another is 
needed to evolve the bra. Those fields are just integrations variables in a path integral. 
Labeling them collectively as 1 and 2 , the full path integral is written as 

In the semi-classical regime the spread of the wavefunction should be small, heuris- 
tically ((f) 1 — <f) 2 ) ~ h((f) 1 +<f) 2 ). We expect this to be connected with the ordinary per- 
turbative regime. To exploit this we make a linear change of variables to the so-called 
Keldysh basis consisting of averaged and differenced fields: 

0^i(0 1 + 2 ), ^ = ^-0 2 . (3.17) 

(We follow the notations used by the finite temperature field theory community.) Now 
we expect simply a ~ h. 

The first step is to express the response functions in this basis. As simple compu- 
tation shows [10], 7 

R n (x 1 ,x 2 ,...,x n ) = -i(O r (x 1 )O a (x 2 )---O a (x n )) (3.18) 

if 0(x\) is the advanced leg. Notice that the perturbations f'O have become the 
differenced operators O a , which happens because it enters with opposite signs on the 
two branches of the contour. 

The next step is to expand the action in powers of <p a . This contains only odd 
powers, the full action S[<p l ] — S[<f) 2 ] being odd under interchange of the two fields. In 



7 When writing O a for a composite operator, the combination 1 ~0 2 should be understood. Thus, 
for instance, for O = (4>) 2 , we let O a = ((j) 1 ) 2 — {4> 2 ) 2 which should not be confused with (4> a ) 2 . The 
same goes for O r = \{0 1 +0 2 ). 
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Figure 4: Examples of tree and one-loop Schwinger-Keldysh diagrams for response 
functions. 

the leading, linear approximation, this becomes a (^f [<p r ] — /) and the (fi a integration 
can be done exactly: it enforces the Euler-Lagrange equations 



Zty° e **°($f[* r ]-/) 



5 



5S 
8<fi(x) 



1-/ 



(3.19) 



Hence, at the leading order, response functions are computed by solving classical 
equations of motion with sources and retarded boundary conditions, as expected. What 
are the corrections to (3.19)? They arise from neglected terms in the full action, 
which are of order (0 a ) 3 . Using (3.18), their effect on some R n is seen to correspond 
to a response function R n +3 with three legs at the same space-time point, where an 
additional interaction vertex is inserted. This has the topology of a 2-loop diagram. 
Hence, the localization (3.19) remains valid at 1-loop! This is a remarkable and far- 
reaching statement, yet it was to be fully expected from the estimate a ~ h. 

Equation (3.19) being valid at 1-loop means that 1-loop corrections only arise 
through initial conditions. The path integral reduces to a sum over classical solu- 
tions, each of which is determined by its own initial conditions, but there may be a 
distribution of such initial conditions. Thus our problem is one in classical statistical 
field theory, and, indeed, (3.19) looks very similar to the Siggia-Martin-Rose formula- 
tion [19] of statistical field theory. At 1-loop it is sufficient to know only the 2-point 
function of the distribution, which is ~ h. The relevant 2-point function is that of the 
averaged field <p r in vacuum, which may be readily computed. For a scalar field 



G rr (q) = n5(q 2 + m 2 ). 



(3.20) 



Inserting this 2-point function into an amputated response function with n+2 legs will 
give the 1-loop response function. For a theory with a real spinless scalar, with proper 
symmetry factors this procedure immediately gives 



R 



! - |o °P( Pl , . . . ,p n ) = - J 7tS(q 2 +m 2 h ) R^ 2 ( Pl , ...,p n ;q, -q). (3.21) 



Putting back spin indices, this precisely reproduces (1.2), providing a very simple and 
physically transparent proof of that result. 
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3.4 Perturbation theory 

Feynman rules in the Keldysh basis are summarized in Appendix A. They are not too 
complicated. For concreteness, here we give the propagator which is a 2 x 2 matrix, 

Q rr Qra \ / I (G > + Q<) G fl \ / nS( P 2 + m 2 ) ^ 

G" G"J \G A )-{^ 

The key feature of this propagator is that G aa = (which is a non-perturbative 
statement known as a "largest-time equation": <p a — y at large times.) Furthermore, 
all vertices contain odd numbers of a fields, and the most "common" vertices have 
only one a field. In practice, this means that almost all propagators in a retarded 
amplitude will be retarded propagators. Heuristically, vertices with more than one a 
are suppressed by powers of h . 

These heuristic estimates can be justified by rigorous counting. Let P denote the 
number of propagators in a diagram, let V be the number of vertices and let n a be 
the number of extra a fields at the internal vertices (n a is twice the number of h - 
vertices). Then, balancing the number of a fields between propagators, vertices and 
external legs (taking n — 1 external legs to be of a- type for a response function), gives 
P — P rr = n a + V + n — 1, where P rr is the number of rr propagators. Combining this 
with Euler's formula for the loop number, £ = P — V + 1 — n, yields 

£ = P„ + n a . (3.22) 

This is the desired formula. Both terms are non-negative and n a is even. Equivalently, 
G rr ~ % and <p a ~ h. s 

Thanks to (3.22) it is a simple exercise to use the Feynman rules of Appendix A 
to reproduce the tree- and 1-loop results obtained above. We leave this exercise to the 
reader. At tree-level, one should find a sum over ordinary Feynman graphs but with all 
propagators being retarded, since (3.22) says that there cannot be any cut propagator. 
At 1-loop, one should get a sum over diagrams that have exactly one cut propagator, as 
for instance in Fig. 4 (b), whose combinatorics are trivially shown to reproduce (1.2). 

4 Formal manipulations of higher-loop integrals 

In this section we derive the results announced in section (2). We begin by systemati- 
cally analyzing 2-loop Schwinger-Keldysh diagrams, which will lead to the formula (2.4) 
involving in general an unphysical sum over tree diagrams, and to its 3-loop extension. 

4.1 Cutting 2-loop diagrams 

The new ingredient at 2-loop is the /i 2 -vertex, that is the vertices with 3 <f) a fields. 
Physically, this vertex is required for unitarity of the theory. Without any loss in 
generality we can simplify the combinatorics by restricting to theories with interaction 
vertices of at most cubic order: Higher-order vertices could be broken down to cubic 



3 For a correlation function involving n r external r- fields, the left-hand side becomes (i + n r — 1). 
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Figure 5: Anatomy of a 2-loop diagram containing a ^-suppressed vertex, shown as 
the dotted one. The retarded propagators on the figure are really strings of retarded 
propagators with insertions of external fields, as illustrated on the right. 

ones at the cost of introducing auxiliary fields; in the Schwinger-Keldysh formalism, 
the propagator of an auxiliary field with kinetic term S = |£ 2 is simply 

G rr G ra \ / % \ ,. OON 

The most general 2-loop diagram with an fi 2 - vertex has the form depicted in Fig. 5. 
The loops need to overlap, with is why this is no other possibility. To simplify the 
discussion we omit the external a fields, which do not play any active role in this story 
- their only effect is to source a classical background field and "dress" the retarded 
propagators in the diagram. 

In terms of dressed propagators, all diagrams reduce to Fig. 5. Let us dissect it 
in some detail. We call the % -vertex "vertex 1". Starting from it and following the 
arrows there are three paths. All three eventually merge together, but two of them will 
always merge beforehand. We label these two shorter segments A and B, ending them 
at vertex 2 where they meet. Then segment C is the third other one, going to the final 
point, and segment D is the merger of segments A and B up to that same final point. 
This labeling, shown in Fig. 5, is unique up to the interchange of paths A and B. 

The obvious symmetry between vertices 1 and 2 suggests a special role for vertex 2. 
Indeed, topologically speaking, in a 2-loop diagram with overlapping loops, there are 
exactly two vertices which can be turned into fr 2 -vertices. These are those two for which 
all incoming momenta are free, not restricted by the diagram's external momenta, and 
are vertex 1 and vertex 2. In that sense they form a pair. 

The propagators in the diagrams are all retarded. To cut-open two of them we use 
causality as explained in the introduction. Here the relevant identity stems from the 
definition of the retarded propagator, 

GR, A (x-y) = ±p(x-y)9(±(x°-y )), 

where the spectral density p(x—y) = [<f>(x) , (f>(y)] is supported, in momentum space, 
on the mass shell q 2 + m^ = 0. Whenever a ^-function of time is redundant we drop 
it and gain a cut propagator. This can be done at the level of "dressed" propagators, 
e.g. strings of bare propagators, using the telescopic identity 

G R G R - ■ -G R = P G R - ■ -G R + G A p- ..G r + ... + G a G a - ■ -p, (4.24) 
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which holds modulo a term GaGa • • • Ga which vanishes when multiplied with a 9- 
function of the time separation between the string endpoints. In words, each propagator 
in a cut string gets cut in turn. Thus we can work with dressed propagators just like 
with ordinary propagators. 

To apply this to the diagram, for instance, one can first cut-open propagators B, 
then use the remaining ^-functions along segments C and (A+D) to cut-open a second 
segment. This leaves a connected tree, and then there is no way to proceed further 
using (4.24). 

We shall need the most general identity obtained by this method and so let us intro- 
duce a hyper-condensed notation. Let G(A) denote the string of retarded propagators 
of segment A and let p(A) is the right-hand side of (4.24) applied to it. (We omit 
subscripts G R) a from the retarded propagators whenever the time flow is clear from 
the context - time always flow from a cut propagator to its future.) The most general 
identity consistent with symmetry between A and B then reads 

G R (A)G R (B)G R (C)G R (D) = a G(A)p(B)p(C)G(D) 

+ /3 p(B)G(C) [p(A)G(D) + p(D)G(A)\ 
+ (A+>B), 

where a = | — (3 is one free parameter. We treat the diagrams where 2 is the h vertex 
in the same way, with the same a, so as to ensure symmetry under 1 -H- 2. Hence 
we have only one free parameter consistent with the symmetries. On the other hand, 
two topologically distinct cuts appear (A-B or A-C), and unique choice will make their 
coefficients equal. This is our best hope at obtaining a "natural" formula. This choice 
is a = jtj and {3 — ^ leading to 

G R (A)G R (B)G R (C)G R (D) + G A (A)G A (B)G R (C)G R (D) 

=i [p(A)G(B)p(C)G(D) + G(A)p(B)p(C)G(D) 

-p{A)G{B)G{C)p{D) - G{A)p{B)G{C)p{D) 
+p(A)p(B)G(C)G(D)} . 

The pattern of signs is the same as in (2.4)! Using the same notation as there, and 
recalling factors of | from the vertices, we have found that the sum of diagrams with 
a h vertex is 



Ii 



2-loop,7i vertex 



i=l 



--- 1 hi J ^ 



E 



^h\ „—hi J12 „—h,2\ 



connected trees 

{1/3, q\ and q 2 parallel, 
— 1/3, qx and q 2 anti-parallel, (4-25) 

0, qi and q 2 non-comparable. 

Integrating out the auxiliary fields to restore original four-point or higher vertices is 
easily shown to pose no problem in this form, but this seems to be the case only for 
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a = j^. Upon including the diagrams with no Ti 2 vertices, which are described by the 



classical contribution (2.3), this establishes the 2-loop formula (2.4). 



4.2 3-loop formula 

The 3-loop case is very easy to analyze once one realizes that a 3-loop diagram contains 
at most one h -vertex. Those diagrams will contain one rr-propagator, which acts like 
a classical fluctuation. Stripping it off leaves a 2-loop amplitude, which we have just 
analyzed. More generally, the same analysis applies to any £-loop diagram with only 
one h 2 vertex. One needs only be careful with symmetry factors. For the sum of all 
3-loop diagrams including those with no U 2 vertex we obtain 



Ti5 + (qi+ml) 



Rl-^( Pl ,..., Pn ) = ijl EH 2 " 1 

i=l L hi J 1 

x Yl R n+e(pi,---,p n ;qi 1 ,- 

connected trees 

2, q\ and q2 parallel, 
x < 0, q\ and g2 anti-parallel, 

1, q\ and q^ non-comparable 



-^ 



■,Qs 3 ,-Qs hs ) 



(4.26) 



At 4-loop level one will find diagrams that contain two % -vertices. Their analysis 
will not be attempted here. 

4.3 Planar gauge theories 

In section 2 we gave a formula for 2-loop corrections in planar gauge theories, which 
we reproduce here 



Ii 



2-loop, planar 



(Pl,---,Pn) 



*m 



i=l 



£< 



\2hi 



n5 + (q 2 +m 2 hi ) 



Y Ru+a{- --,Pi, gi, -qi, ■ ■ ■ ,Pj, 92, -92, • • •) 

.i<3 

+ X^ ^™+4(- • ' >P«> 9l, -92, 92, -9i, • • •) 
i 

4/3, sgn $ = sgn q%, 
2/3, sgn g° = -sgn q%. 



(4.27) 



We now show that this Ansatz follows from (2.4), which is true in any theory. 

Indeed, (4.27) and (2.4) only differ in the treatment of "non-comparable" Feynman 
diagrams, which are now made "comparable" by their color structures. We have to 
show that this difference integrates to zero. 

The relevant cancellation pattern is embodied in (4.24). It is a multi-particle gen- 
eralization of (4.28). Consider first the disjoint color structure, the first case in the 
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Ansatz. In the "non-comparable" diagrams, by definition, the q± and q 2 loops are dis- 
joint from each other. (They may share a vertex, which will not interfere with the 
present argument, but they share no propagator.) Each such Feynman diagrams will 
generate one term for each choice, in each loop, of which propagator is cut. The color 
flow and the momentum flow follow each other along the loop so the sum is of the 
telescopic form (4.24). But now we have a telescopic sum around a full circle, so both 
qi and q 2 integrands reduce to products of retarded propagators. They integrate to 
zero. 

The "nested" color structure is treated similarly. The Feynman diagrams in this 
class that have two disjoint loops now simply have tadpole insertions; this class of 
diagrams does not overlap with the ones just considered. 

Equation (4.27) is still a rather formal result. Difficulties appear in the forward 
limit q' — > q: Tadpoles divergences 

{q - q'Y 

no longer cancel, as the familiar cancellation in gauge theory requires a plus sign 
between the two terms. In field theory parlance, the familiar cancellation of ~ f abc 
tadpoles when contracted with 5 ab is spoiled by the combination (2.7). This means 
that in planar gauge theories one trades the unphysical amplitude (2.4) to the forward 
limit of a physical one, at the cost of a dependence on how the forward limit is taken. 
The proper analysis of this ambiguity lies beyond the scope of this paper. We mention, 
however, that for sufficiently supersymmetric planar theories the tadpoles cancel in 
the susy sum and (4.27) is unambiguous as it stands! This will be discussed further in 
section 6. At the level of an individual diagram, the equation is tested in Appendix D. 

4.4 Evidence for the uniqueness conjecture 

Finally we provide evidence, actually sketch of a proof, starting from (4.25) of the 2- 
loop uniqueness conjecture (2.10). The setup now is a general non-planar quantum field 
theory. We will argue that the integral of the quantum term (4.25) is fully determined 
by 3-particle unitarity cuts, only, evaluated at complex loop momenta. We will then 
argue that these cuts are physical and determined by their classical counterparts, which 
together would imply that a loop integral with vanishing classical 2-particle cuts is zero. 

The strategy is to employ contour integration in the rapidity plane of one of the on- 
shell particles. This requires space-time dimension D > 2. When D > 2, the choice of 
a rapidity plane breaks Lorentz invariance, and the usefulness of the resulting formula 
for practical computations is unclear. It suits, nonetheless, our present purpose. 

The cancellation which will play a role is exemplified by the following integral: 

[ n6(q 2 )sgn(q°) [G R (q+p) + G A (q-p)} = 0. (4.28) 

Jq 

This integral vanishes by construction because, modulo shifts of q, it equals the contour 
integral of ^-/q 2 {q+p) 2 at infinity in the complex frequency plane. The only reason the 
integrand in (4.28) appears to be non-zero is the relative shift between the two terms. 
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Figure 6: Integration contour C in the rapidity plane used in sectin 4.4. The interior 
of this contour is characterized by Imp z > 0. 

Really, the integrand is zero, this is just hard to see. More precisely, the bracket 
times the sign function is zero when viewed as a function on phase space and a certain 
relative shift between the two terms is made. In any case, this vanishing will become 
more manifest once we perform a second contour integration, which will break the shift 
symmetry. 

To do that second integration, the crucial point is to note that (4.28) is itself a 
contour integral in the phase space parametrized by q. This crucially relies on the sign 
function. To see this, introduce a rapidity variable such that g M 
^ q^_+m 2 sinh r] , q±). The integration measure becomes 



( ^/qj_+m 2 cosh 77, 



/ ' 7v5(q 2 )sgn(q°) 
J a 



d g_L <p drj. 



The contour C has two real segments located at \mr] = 0, ir, as depicted in Fig 6, 
corresponding to positive and negative energies. Crucially, the sign function provides 
the correct orientation for the contours. The contour can be closed at infinity on each 
side by adding small vertical segments, the integral over which vanishes. Hence this 
integral is determined by the residues in its interior, which characterized by Imp 2 > 0. 
(This "interior" may not appear parity-invariant. However, the lines Imp z = 0, on 
which the original contour runs, divide the complex phase-space into two halves which 
are parity conjugate to each other. Closing the contour on either half will give the 
same answer, which is thus parity-invariant.) 

If the first term in (4.28) has a pole at q — q, then the second term has a pole at 
q — p corresponding to the same physical configuration, e.g., the same unitarity cut of 
a bubble integral. If we assume that Imp 2 = (which we certainly can, for generic, 
off-shell external momenta), the two poles are always inside or outside the contour at 
the same time. The residue of the first pole is 



pOqz _ pZgO 



The residue of the pole in the second term is 

-1 



p°(q z —p z )—p z (q°—p ) 
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Since the shift does not change the answer, the residues are equal and opposite and we 
see that (4.28) vanishes, as expected. 

This cancellation pattern, for phase space integrals involving sgn(g°), is all we need. 
We apply the same technique to the qi integral in (4.25). Let qi be a cut propagator 
along, say, segment B (in the notation of Fig. (5)). The two extra propagators cut 
by the q\ contour integral can then be on AA, CC, DD, CD, AC or AD. However, 
all but the last two cases AC and AD are easily seen to cancel by the above pattern. 
The remaining cases AC and AD correspond, including cut propagator along B, to 
genuine 3-particle cuts, e.g. the exchange of a 3-particles state, and the reason they 
do not cancel is due to a relative sign introduced by the weighting in (4.25). The same 
conclusion is obtained regardless of where q2 is in the diagram: only genuine 3-particle 
cuts survive. 

There remains the claim that the required 3-particle cuts are just ordinary uni- 
tarity cuts and so are determined by (residues of) the forward limit of physical tree 
amplitudes. It is hard to see how this could not be the case, given that in the real 
but unphysical region considered around (2.5), perfect cancellation must occur between 
the amplitude needed here and a (residue of a) physical forward amplitude. Then by 
analytic continuation this should continue to hold at complex momenta. It would be 
good to see this more explicitly. 

5 Example: 2-loop vacuum polarization 

To help make these considerations more concrete, let us now compute the 2-loop vac- 
uum polarization in massless QED 

Il£» = -i I d 4 xd- ip - x {0\[r(x),J u (0)}\0}6(x°). (5.29) 

The current will be off-shell, p 2 ^ 0. This quantity was computed long ago by Kallen 
and Sabry [20]; the massless limit enters the the 3-loop electron contribution to the 
muon g— 2. It also appears in QCD, where its imaginary part is related to the leading 
0(a s ) corrections to the total e + e~ — > hadron cross-section. We will see here that the 
well-known result can be derived using only on-shell data, even though we make no 
claim as to the efficiency of this approach for this particular computation. 

5.1 1-loop vacuum polarization 

Let us begin at 1-loop. We need the forward amplitude of a quark against two currents, 
as in Fig. 7 (a). Labeling the momentum of the incoming quark as q the forward 
amplitude is 

A^(q) = e 2 Tr [w^^-f + ^(^^] ' ^ 

where q 2 = 0. The helicity sum in (1.2) simply adds the same thing with q — > —q. We 
have absorbed the minus sign (— ) 2h into (5.30), so (1.2) becomes 



n£T(p) = /WMfte). (5.3i) 

J Q 



23 



■X- 



-X — X 

(a) 



X X 



■X- 



(b) 



(c) 



Figure 7: Forward amplitudes against two off-shell currents, denoted as crosses, con- 
tributing to 1- and 2-loop vacuum polarization. All external legs are on-shell. 



The integrand is manifestly transverse: p^Ai U (q) = 0. Transverseness tradition- 
ally holds only for integrals, but, in an on-shell formulation, it is already manifest at 
the integrand level. Lorentz-covariance is also manifest. Transverseness and Lorentz 
invariance imply, in the standard fashion, that we only need to compute the trace. 
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(5.32) 



In the bracket we have dropped a tadpole not depending on p, which vanishes in 
dimensional regularization. 

This D-dimensional integral may be readily carried out like the previous triangle 
example. An equivalent, but much simpler, way to proceed is to recognize (5.32) as 
the on-shell representation of a scalar bubble, and evaluate the latter using standard 
Euclidean space techniques. That is, 
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(5.33) 



The first equality holds because taking single-cuts of the right-hand side yields the left- 
hand side. (In Minkowski signature for the final result one has to remember the branch 
of the function (p 2 )~ e appropriate for retarded boundary conditions, p 2 — > p 2 
This way the well-known MS result is reproduced: 
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5.2 2-loop vacuum polarization 

In the present example a new subtlety appears at 2-loop: pinching singularities oc 
S(q 2 )/q 2 due to self-energy subdiagrams. Indeed the formula (2.4) contains no instruc- 
tion to "amputate," and, besides, it is hard to see how one could "amputate" diagrams 
in the present framework without spoiling gauge invariance. 
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We postpone a more satisfactory treatment of these singularities to subsection 5.3. 
For the time being we just use a quick fix. Thinking of the function 5(q 2 ) as the result 
of a contour integration over a small circle around 1/q 2 , suggests viewing S(q 2 )/q 2 as 
a poor guess for a contour integral around l/2g 4 . Computing a residue around l/2g 4 
ought to be the right thing to do. This leads to the prescription, for any F(q) regular 



at q 2 = 0, 



^F (q) - I 
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F{q) 



2(q 2 - ie) 



-n5'(q 2 )F(q). 



(5.35) 



The symmetry factor 1/2 is because when two propagators pinch in a diagram, S(q 2 )/q 2 
gets generated once for each propagator. Treating each propagator using (5.35) will 
produce the right answer. Using this prescription we can evaluate the forward ampli- 
tude of a fermion against two currents plus two photon insertions, depicted in Fig. 7 
(b) (summed over permutations), which contribute to II^ an amount 



- J \s'(q 2 )7r8( q 2 2 )Tr 



1 



1 



9i e 



9i + 92 



e 9i7" 



1 



9i +P 



Y ' + ((p,ij)++ i-p,")) 



l 



9ie 



-7 



7 



9i + 92 9i + 92 + P 9i + 92 



+917' 



+qiY 



-r 



9i + P 9i + 92 + P 9i + 92 

1 1 1 

e e* 7* 

9i + P 9i + 92 + P qi+P 
1 ,. 1 1 



+ qi e Y = e*^— Y 

9i + 92 9i + 92 + P Qi+P 

+(G»,/i) ++(-p,v))]. 



(5.36) 



This is to be summed over photon polarizations e; we will systematically use symmetry 
under q — > —q to simplify the writing of the integrands. Keeping this in mind and using 
the identity q 2 5'(qf) = —S(qf), it is easy to check that (5.36) is manifestly transverse 
under both e a — > q^ and 7 M — > j>. 

It is important to stress that (5.36) includes only physical polarizations e, e* for 
the exchanged photon. No gauge fixing is needed because that photon is on-shell. 
Longitudinal photon polarizations manifestly decouple at the level of the integrand. 
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This implies that the polarization sum is equivalent to a Lorentz trace ^2 

This fact will be used here, for convenience, as it helps simplify the evaluation of the 

individual terms as opposed to their sums. 

The trace of (5.36) is still a complicated function. To avoid unimportant details, 
we now show details only for the terms having one particular Dirac structure, Tr [-/ty M -^ 
lti(py a p~Ya], with obvious notation, corresponding to the first, second and fourth line 
of (5.36). The mass-shell delta-functions allow the numerators to be simplified and 
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II R ^/e 4 to be reduced to 
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Next we need to compute the forward amplitudes of two fermions against the cur- 
rents as depicted in Fig. 7 (c), which are listed in the Appendix C. These are manifestly 
transverse by themselves and again it suffices to compute their trace. The terms that 
have the same Dirac structure as above are 
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(5.31 

The cancellation of the 8' terms between (5.37) and (5.38) is due to chiral symmetry, 
which forbids a massless electron from acquiring a mass. 

At this stage it would seem that we have many non-standard integrals to com- 
pute. Furthermore, we don't even have the whole integrand yet since the non-standard 
amplitudes with the |'s in (2.4) are not yet included. 

The beauty of the formalism is that this is unnecessary. The on-shell integrals 
are guaranteed to sum up to the classical cuts of standard Feynman integrals. In the 
present case, it is easy to see that modulo shifts in the integration variables, the sum 
of (5.37) and (5.38) is reproduced by the classical cuts of 
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(5.39) 



Xqi+q2) 2 (qi+p) 2 (qi+q2+p) 2 . 

This is now a set of standard Feynman integrals. 

It seems intuitively clear that the passage from (5.37) and (5.38) to (5.39) is unique, 
modulo equivalent ways of rewriting the latter. In general, this is the content of the 
conjecture (2.10). In the present case there is actually nothing conjectural about it, 
since the exact same computation could have been done in a planar gauge theory. 

Repeating the same steps for the remaining Dirac structures in (5.36) and (C.60), 
the forward amplitudes are found to uniquely determine the vacuum polarization as 
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(5.40) 
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The first term may be readily integrated using integration-by-parts techniques of [21] 
and the other integrals are straightforward. 

Actually, the first line is just what one would have obtained from computing ordi- 
nary Feynman diagrams, as it must. The point of the present exercise was to deduce 
it directly from its cuts. 

As an extra consistency check on the formalism, the scalar integral corresponding 
to the second term in the bracket is evaluated in its full on-shell form in Appendix (D), 
giving a nontrivial check of the factor | in (2.4). The result agrees. 

The imaginary part of (5.40) gives, with the proper color factors, the famous cor- 
rection [1 + ^-Cp] to the total e + e~ hadronic cross-section in QCD. We deduced it 
here starting from physically observable forward scattering amplitudes. 

5.3 Pinching poles and S'(q 2 ) from forward limits 

To complete our claim that (5.40) comes from physical data, it is important to check 
that the S'(q 2 ) terms in (5.36) does come from physical data. This is not entirely 
obvious, because q 2 S'(q 2 ) ^ would suggest that the formula actually needs informa- 
tion about infinitesimally off-shell particles. We show here that the required data is 
contained within the forward limit of a physical on-shell amplitude. 

We consider the following on-shell, but non-forward, deformation of the tree ampli- 
tude: 

-Rn+4(e) = Rn+i{P - e<?2, ~P\ <?1, -<?1, <?2 + ^2, -^)- 

This conserves total momentum. The limit e — > has to be taken carefully for the 
Feynman diagrams of self-energy type, for instance Fig.l (d) with p -f-> q 2 . In such 
diagrams, the two terms in which the cut fermion is on either side of the self-energy 
insertion sum up to a 0/0 form 
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Here Al and Ar stand for the two "self-energy" factors on the two sides of the singular 
propagator. Now using the identity 

ri(q>6(q 2 2 )-^J^F( qi , q2 ) = - [ xfi&KSi&Ffato), 
iq\-q 2 oq 1 



qi,Q2 "'^ i 'i^ "11 ^ 91,92 



we see that the limit e — > exactly reproduces the first line in (5.36). This establishes 
that only physical data enters (5.36). 9 



'Sometimes it may be useful to have a deformation preserving a mass-shell condition for all legs, 
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6 Supersymmetry 

Supersymmetric theories are special. On a somewhat more modest scale, this is also 
seen when the present formalism is applied to them. 

6.1 LSZ reduction 

The above-described formalism has been set up for the computation of off-shell cor- 
relation functions. Ultimately, we want to consider on-shell S-matrix elements, which 
is where its full power should become manifest. To do so, the LSZ reduction pro- 
cedure has to be adapted to make it applicable to loop integrands. As a first step 
toward this goal we now present how this works for sufficiently supersymmetric theo- 
ries, where supersymmetry cancellations simplify this task. Here "sufficient" will turn 
out to mean M = 1 for massless particles and M = 2 for massive particles. Discussion 
of non-supersymmetric theories is left to future work. 

To see the subtleties, consider an external gaugino leg, in M = 1 supersymmetry, 
whose momentum p approaches its mass shell p 2 = 0. In this limit the forward tree 
amplitude contains terms in which in which the loop momenta q and — q are inserted 
on the p leg, as in a self-energy insertion, leading to terms which diverge like 1/p 2 . 
Normally one would amputate this a-la LSZ, defining 

S( Pt ) = (Tl t Zl /2 p 2 )A( Pt ), 

where the Z$ are appropriate wavefunction renormalization factors, S(pi) is a S-matrix 
element and A(pi) is the sum of graphs for a correlation function omitting self-energy 
subgraphs. Now in the present context it is not clear how to "amputate," especially 
if the data one is given consists of the physical forward amplitudes summed over all 
graphs, as we assume. 

The singular terms arise from well-separated self-energy subgraphs: thus the anal- 
ysis boils down to a 4-point function. Including only the neighboring propagators, the 
forward amplitude of a gauge boson against our gaugino reads 



A 



i>+ i j* + j* p 1 - 



(p+q) 2 ' (p-q)' 
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where (e\,q) describe the states of the (forward) gauge boson. 

Upon summing over gauge boson polarizations one finds the divergent contribution 
to the forward amplitude 



-2g 2 C A x 1 (6.42) 

^4 



for instance p 2 = 0. In this case one could consider a more general deformation such as 

R €l ,€ 2 = R(jP + eitfi + £ 2<?2, -p, qi + eigi, -qi,qz + e 2 q 2 , -qi-, ■ ■ ■) (5-41) 

where e\ = —(p-q2)/(p-qi)^2 and the dots stand for an arbitrary number of other external states. We 
find that this can also be related to (5.35), upon subtraction of certain gauge-invariant terms involving 
the scale dependence qi^ir of the (physical) on-shell self-energies. 
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Such a divergence is interpreted as a gaugino mass being generated. This may 
seem puzzling, given that this should be forbidden by gauge invariance and Lorentz 
invariance. However, one cannot invoke Lorentz invariance at this stage, since the q 
integral is not performed yet! In fact, there is no way the forward amplitude could 
have vanished at the integrand level: in the presence of a thermal gas of gauge bosons, 
the gaugino would rightfully acquire a (non-Lorentz-invariant) thermal mass, and it is 
precisely the job of the forward amplitude to generate it [22]. What happens in vacuum 
is simply that (6.42) integrates to zero by virtue of Lorentz invariance. (The integrand 
is a total derivative, J S(q 2 ) oc J -^q tl 5{q 2 ) = 0, and, as such, it should vanish in any 
reasonable regularization.) 

Gauge invariance does not suffice to obtain cancellations at the integrand level. 
Supersymmetry, however, provides such cancellations. In our case, the superpartner of 
the gauge boson produces the forward amplitude of a gaugino described by (X q , q) 
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This diverges like 
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This exactly cancels the bosonic divergence, so that the sum possesses the milder 1/p 2 
behavior required for the LSZ amputation procedure. 

It can be seen that these cancellations are general features of supersymmetric the- 
ories. The point is that divergences are proportional to 1/p 4 times physical, on-shell 
gauge-invariant four-point scattering amplitudes. Such forward four-point amplitudes 
are just c-numbers, since all non- vanishing Mandelstam invariants are equal. In super- 
symmetric theories one has to sum over the polarizations of two of the four particles, 
q and —q. The sum will vanish provided the other two states, p and —p, form a BPS 
object. For a massless particle in M = 1 this is the case (the two supersymmetry 
transformations X"Q a and X^Qa annihilate both states), and for any particle in M = 2 
supersymmetry this is also the case. Equality of asymptotic thermal masses in a finite- 
temperature context for different particle species occurs by a similar mechanism [23]. 
In all cases here, the mass shifts cancel. 

After the leading divergence is so canceled, it is important that the subleading 
term be a clean 1/p 2 , as required for LSZ reduction. This is easily verified to be the 
case quite generally, as a consequence of the Lorentz-invariance of the equations of 
motion of the theory - the self-energy sandwiched between propagators (which provide 
j> polarization tensors) has to vanish when p 2 = 0, and after 1/p 4 is canceled we find 
that one is always left with a clean 1/p 2 . 

In summary, in sufficiently supersymmetric theories the following limit exists 
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The wavefunction factors Zi can be defined through a similar limit 
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where G4 is a four-point correlation function. The combination of the two objects 
produces correct S-matrix elements S (with retarded boundary condition). 

The gauge dependence in Zi and in the on-shell limit of R cancel against each other, 
making S gauge-invariant as usual. In fact we have checked, using Ward identities, 
that this cancellation does takes place already at the level of the forward amplitude - 
longitudinal gauge bosons and ghosts cancel and need not be included in the helicity 
sum. 

6.1.1 Forward limits 

The following method might be a convenient alternative way to evaluate the forward 
amplitudes. The reader should be aware, however, that we have only tested it in M = 4. 
This method takes the external momenta to be on-shell from the start, which makes 
gauge invariance fully manifest, but takes the amplitude to be slightly non- forward 

lim J2(-) 2hRt ™2(Pn •••»»,, q\ (-q+e^)- h ) (6.45) 

h 

where q-k = and now R n +2 is already amputated. It is easy to deform the other pi 
so as to conserve total momentum and simultaneously preserve mass-shell conditions 
pf = —m 2 , as in section 5.3; because mass shift divergences now cancel any such 
choice is equivalent. A-priori, now, this limit may depend on a arbitrary relative phase 
between the polarization tensors of q and (—q+ek). This will disappear if complex- 
conjugate polarizations are used for CPT-conjugate states, as is natural: then the 
ambiguity will be in the cosine of a phase and drop out. Provided the /x-dependence 
drops out, which we have explicitly verified in M = 4, this gives a reasonable candidate 
for the correct answer. 

6.2 Good ultraviolet behavior and factorization from unitarity cuts 

Consider deforming a forward amplitude in a gauge theory (we will consider only non- 
gravitational theories, for simplicity) by taking q\ — > q± + z\x where qy\x = fi 2 = 0. 

R n+2 (z) = Y^(-) 2h R n+ 2(pi, ■■■,Pn,(q + zn) h , {-q - zfi) h ). (6.46) 

h 

Here again R n +2 is already amputated and all momenta are on-shell. This defines a 
rational function of complex z whose poles are due to internal propagators becoming 
on-shell. This is similar to a Britto-Cachazo-Feng-Witten (BCFW) deformation [24], 
but somewhat simpler because really only one momentum is shifted. As in the BCFW 
case, R n+ 2{z) is completely determined by its poles provided it vanishes at infinity, 
which leads us to consider its behavior at infinity. 

Direct inspection of individual Feynman graphs reveals that R n+ 2{z) grows at most 
as z at large z. This is the familiar eikonal limit. (Here polarization tensors can 
be taken to scale like ~ 1 in this limit.) However, since the eikonal contribution is 
spin-independent it cancels among superpartners, making R n+2 (z) is at most of order 
one. 
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To find the true behavior of R n+2 (z) a more refined analysis is needed. Supersym- 
metry cancellations can be enforced from the start by using M = 1 supersymmetry 
to trade the sum over superpartners to a sum over supersymmetry-transformed back- 
grounds, 

A O a i 

h bosonic 9 

where Q acts only on the p states; A is an arbitrary reference spinor such that the 
denominator is nonzero, where X q (z) is defined by $(z)\ q (z) = 0. This is now a 
"forward" amplitude in which a boson turns into a fermion. The leading diagrams 
have one gaugino exchange and a certain number of gauge boson exchanges. (F-term 
interactions are not relevant because they couple, for instance, a holomorphic scalar 
to holomorphic fermions ip, while it is the complex conjugate fermion ip that appears 
in the forward amplitude.) In fact, gauge bosons can also be neglected since the light- 
cone gauge q-A = is non-singular for the relevant diagrams. Hence the large z limit 
is contained in a single tadpole graph involving the exchange of a zero-momentum 
gaugino, and proportional to Trt a . 

Such a tadpole graph, if nonvanishing, would of course be ill-defined. In fact, when 
such tadpoles graphs with zero-momentum propagators exist, the forward amplitudes 
themselves are ambiguous to begin with, and what we are seeing at large z is just 
this ambiguity. This possibility was ignored in the last subsection; in supersymmetric 
theories it can happen only if there is a gravitational anomaly (Trt a ^ 0). Ignoring 
this possibility and assuming the absence of a gravitational anomaly, then forward 
amplitudes exist and the above argument shows that they vanish at infinity. 

This implies that forward amplitudes can be constructed from their poles: 



Rn+2(q) = Yl / dr h dr h a l({q + ^(i>, m}, {pi -(q + z ( I )v),m}, m L ) 
i •* 

;M{-(<1 + z(I)»), -Vi}, i-Pi + (Q + <!)»), -V2}, Vi R )- (6.47) 



(pi+qf 

Each term is associated to a 2-particle factorization channel I with total momentum pi, 
that is a 1-loop unitarity cut, with z(I) = — Vl 2 q Pl corresponding to the point where 
(pi + q + zfi) 2 = 0; the 77 variables are Grassmann variables to do the sum over internal 
states. This nicely expresses forward amplitudes (in either D or 4 dimensions), and thus 
loops, in terms of the integrands of unitarity cuts (in D or 4 dimensions, respectively). 
No (bosonic) integration over these unitarity cuts is needed to find the integrands. 

Although we did not make use this fact here to obtain (6.47), it may be useful to 
note that in supersymmetric theories only 4-dimensional integrands are needed [2]. 

Two classes of theories have been left apart in the preceding two subsections: theo- 
ries with a gravitational anomaly, and, more interestingly, M = 1 theories with massive 
particles. On the other hand, (6.47) would appear to provide a definition of forward 
amplitudes in these two cases. It would be interesting to further investigate this pos- 
sibility. 
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6.3 J\f = 4 SYM examples 

It is illuminating to illustrate these findings in the most supersymmetric field theory, 
M = 4 SYM. For the present work we shall content ourselves with two simple 1-loop 
examples. 

6.3.1 4-point function 

Consider the 4-point amplitude i? 4 (l + , 2 + , 3~, 4~) in planar M = 4 SYM. Four different 
color structures contribute in the planar limit. Let us focus on the one involving 
£ 6 (-g,g,l + ,2+,3-,4-). 

The factorization (6.47) of this amplitude yields 3 terms. It would be possible to 
evaluate them, but, actually, in M = 4 SYM, a much simpler factorization exists. This 
factorization involves taking \x — q, e.g. deforming q — > zq. In general this would lead 
to double poles at z = (corresponding to the familiar soft divergence ~ l/q-piq-p^), 
which is not very convenient. However, in M = 4 SYM, a simple argument given in 
Appendix E shows that the amplitude decay like 1/z 3 at large z. This is in keeping 
with the finiteness of the theory. 

The "bonus" obtained from this good decay is that one can factor the function 
z 2 Rs(zq) instead. This function has no pole at the origin, leaving only a single channel 

k ^) = i T Z \ ,2 EH^leftMiW^), (6.48) 

{q+Pi+P2) 2 ^ 

where Zq solves (^o9+Pi+P2) 2 — and both factors are MHV four-point functions. 
Recalling the MHV four-point function 

^' 1+ ' 2+ ' r - 4 '»- (12)(23) 3 (34)(41> - (6 - 49) 

(and using its simple supersymmetric counterparts) we have 
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where p 5 is the on-shell momentum z q+p 1 +p 2 . Here (ij) = e^XfXj is the ordinary 
spinor product with e 12 = 1 the antisymmetric symbol. (Null four-momenta are ex- 
panded as p^ a = \f\f-) Using the Schouten identity, the numerator becomes simply 
(23) 4 (g5) 4 . Finally, using that (5| oc [q(pi+P2)\, we obtain 

(23) 4 —st 

Reiq) = (12) (23) (34) (41) X ( q+Pl )2( q+Pl+ p 2 y( q -p 4 y (6 - 51) 

The first factor is the tree amplitude. The three other color structures are computed 
similarly. Summing them, we "discover" that 



,1-loop Atrce . . 2 



^- 10 °P = i^ce x _ g 2 Nst I nS ( g 



+ cyclic 



(q+Pi) 2 (q+Pi+p 2 ) 2 )(q-P4) 
32 



The sum is over the four cyclic permutations of (1234). 

In this expression we recognize the scalar box in disguise: 

«J- = *- x -?S« l qHq+pinq+ ; i+pinq _ n)r (6.52) 

That way, the present formalism quite simply leads to the well-known answer [25, 2]. 
The integral will not be evaluated here, as we keep our focus on the derivation of the 
integrand. Note that (6.52) was obtained here without any a-priori assumption about 
its form. We did not even assume that it could be expanded in terms of boxes. 

A clarification is in order. While we have evaluated the scattering amplitudes 
in 4 dimensions, the integral is .D-dimensional. As explained in [2], neglecting the 
dependence of the integrands on the (— 2e) extra dimensions is justified whenever this 
dependence does not overlap with ultraviolet divergences, which is totally expected in 
M = 4 SYM. 

6.3.2 5-point function 

Let us present a slightly less trivial example, the 5-point MHV amplitude. The inte- 
grand may always be written as 

pl-loop r »/- 

' - S 2 - „_,_»,_ „, „. .„ (6-53) 



R£ ee J q q 2 (q+Pi) 2 (q+Pi+P2) 2 (q-P5-P4) 2 (q-P5) 2 

where M is a polynomial in the loop momentum q. We will now determine M from the 
forward amplitudes. M does not depend on external polarizations, as MHV amplitudes 
with different polarizations are related to each other via supersymmetry. 

Consider the forward amplitude in the color channel between particles 5 and 1, 
summed over a M = 4 multiplet of states. This can be obtained from a 2-term expres- 
sion similar to (6.50) or by any other means, giving 



R 



tree 



(g+Pi2) 2 g54-S5i + (q-Pb) 2 S2iS 2 3 - Tr L [qpmpmPb] 
(q+Pi) 2 (q+Pi2) 2 (q-P45) 2 (q-Pn) 2 



(6.54) 



Here Tr L [qpiP2P3P4P5] = (si) [12] (23) [34] (45) [5?] (in our convention in which the MHV 
tree amplitude depends only on angle brackets). Now this expression is defined only 
when q 2 = but it must agree with the residue of (6.53) at q 2 = for any such q. This 
fixes M up to an ambiguity proportional to q 2 : 

H = (q+Pi2) 2 s 5 4S 51 + (q-p 5 ) 2 s 2 is 2 3 - Tr L [qpmPzPiPs] + q 2 C x . (6.55) 

The ambiguity C\ is ^-independent because the single cuts in other channels must 
decay like 1/q 3 . The same computation in the channel between particles 1 and 2 shows 
that 

N = (q-P45) 2 Si 5 Si2 + q 2 s 3 2Su - Tr L [qp 2 P3P4P5Pi] + (q+Pi) 2 C 2 (6.56) 

where C 2 is another g-independent coefficient. Equations (6.55) and (6.56) must both 
hold for any value of q, and together they fully determine M. Indeed, (6.56) uniquely 
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determines C\ and, vice-versa, (6.55) uniquely determines C%. After a short bit of 
algebra one finds 

N = (q+Pi 2 ) 2 s 5 4S 51 + (q-p 5 ) 2 s 2 is 2 3 + ^q 2 s 32 s 34 ~ ^ L [qp 1 p 2 p 3 p 4 p 5 ] - -g 2 Tr [<y 5 p 2 p 3 p4P 5 } • 

This together with (6.53) constitute our final expression for the 5-point MHV 1-loop 
integrand. By construction, it is rigorously equal to the gauge-invariant sum of order 
~ 100 1-loop Feynman diagrams, modulo purely algebraic manipulations. 
It is possible to rewrite this result in somewhat more symmetrical form, 

N = -g 2 S32S34 + ^(q+Pl) 2 S i3 S45 + ^(q+Pl2) 2 S 5 4S 51 + -(g-J5 45 ) 2 Si 5 Si2 + -(q-p5) 2 S 2 lS 23 

I I 

-2 Tr [l5qPlP2P3P4P5} ~ 2<? 2Tr [75P2P3P4P5]- (6-57) 

The first line gives a sum of five 1-mass boxes and agrees precisely with the form of 
the integrand obtained long ago in [2], using unitarity-based methods. 

The second, parity-odd, line is not manifestly cyclically invariant; it nevertheless 
is. Its overall sign would be opposite for a MHV amplitude. Using Passarino-Veltman 
reduction, however, it is possible to show that it integrates to O(e), which is why this 
contribution is usually not quoted. In that sense, integrands (and their single-cuts) 
contain more information than any of their integrals. 

We have similarly checked that the formalism works for general n-point MHV am- 
plitudes. 10 

It is striking that forward amplitudes in only two channels "know" about the am- 
plitude in all other channels. The magic, evidently, lies in the fast decay of forward 
amplitudes in M = 4 SYM. Using this decay ~ l/q 3 , actually, it seems very likely that 
one could iterate the thinking underlying (6.47) to determine individual the BCFW 
terms in terms of triple-cuts, which would then be determined by quadruple-cuts, at 
which point one would run out of decay properties. Thus forward amplitudes in M = 4 
should be determined by the leading singularities of the theory 11 For this reason we 
expect that, in general, forward amplitudes in only (n— 3) channels should suffice to 
determine n-point 1-loop M = 4 integrands. 

7 Summary and outlook 

Exploiting the physical principles of causality and locality, we have obtained relations 
between loop amplitudes and the forward limit of tree amplitudes, which hold for 
correlation functions in any quantum field theory. At 1-loop the relations are known 
as Feynman's tree theorem, (1.2). Use of causal boundary conditions has allotted us 
significant simplifications to the formulas and discussion, while sacrificing none of the 



10 In that case, in addition to the familiar sum over 2-mass easy boxes, the formalism also produces 
the parity-odd pentagons which were discussed in [3]. Although these pentagons integrate to zero, 
they cannot be dropped through any sequence of algebraic manipulations. 

11 More precisely, the data here would be the leading singularities of the integrand on all possible 
contours, as opposed to only the parity-even combinations which are sometimes considered. 
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physics. A simple analytic continuation relates retarded amplitudes and time-ordered 
amplitudes, but retarded amplitudes could also be useful on their own for computing 
inclusive observables [26]. 

The generalization of the tree theorem to higher-loop is nontrivial. Only for planar 
gauge theories does one gets physical tree amplitudes. In that case our results allow the 
interpretation of the multi-loop vacuum as a collection of on-shell particles with com- 
pletely classical behavior. In the infinite past (where interactions are turned off), these 
particles exist in an intricately entangled but universal color state. We have computed 
this state up to the 3-loop approximation but we attempted no 4-loop computation. 
Already at 2-loop, this color state involves interesting factors of | as given in (2.7). 

In non-planar theories a similar picture holds at 1-loop but not to higher-loop. 
Our discussion of non-planar theories is still rather preliminary and we leave it in a 
somewhat puzzling state. On the one hand we presented a strong argument in section 
2, that integral representations in Minkowski space involving physical forward ampli- 
tudes do not exist then. We believe that this argument is insensitive to the boundary 
conditions used, retarded or time-ordered. On the other hand, we proposed that loop 
integrals are nevertheless implicitly determined by physical forward amplitudes. This 
conjecture, if true, would allow for a purely algebraic determination of higher-loop 
integrands starting from physical forward amplitudes. That would be complimentory 
to unitarity-based approaches, in that one also starts from an Ansatz for the loop 
integrand and then fixes the coefficients. The proper physical interpretation of the 
conjecture, if true, is still unclear. Our support for it comes from a computation in 
complexified Minkowski space-time that we presented in section 4.4. 

The conjecture is also supported by the fact that it would be implied, if unitarity 
indeed suffices to determine loops. Since unitarity-based techniques have been success- 
fully applied to even 4-loop in a non-planar context (Af = 8 supergravity) [27], this 
is actually fairly strong support. It is important to stress that these considerations 
are rather formal. In the end, the only real test of the conjecture will be to work out 
explicit non-planar higher-loop examples. 

We have deliberately left out of the discussion the computation of forward ampli- 
tudes for S-matrix elements. For off-shell correlation functions, to which the formalism 
applies, we believe that everything is under control. But for S-matrix elements, prob- 
ably the most interesting observables, LSZ reduction at the integrand level must still 
be thought through carefully. For supersymmetric theories we did work it out fully, in 
section 6, and this gives us hope that technical difficulties in the general case will be 
solved in the future. 

Nowhere in this paper did we use that the theory was in a Lorentz-invariant 
state. For this reasons, all our results extend to finite-temperature contexts. Since 
we used the Schwinger-Keldysh formalism this is actually quite straightforward - 
formally one just needs to add physical particles to the vacuum ones, specifically 
±7iS(q 2 ) — > 27T<5(g 2 )[±| + n(q)], where n(q) = l/(e q l T =F 1) are particle distributions 
as usual, in (1.2) and in the classical part of (2.4) at 2-loop. In that case our formulas 
are conceptually very close, albeit distinct, from other formulas such as the celebrated 
Dashen-Ma-Bernstein formula for the thermodynamic potential [28] and also the cor- 
relation functions formulas in [29, 30]. These formulas involve vacuum loop amplitude, 
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while the present ones involve vacuum tree amplitudes. 

Finally, it is worth stressing that for sufficiently supersymmetric planar gauge theo- 
ries, no technical subtleties arise and the formulas in this paper directly give loop-level 
S-matrix elements, at least up to 3-loop, as integrals over trees. These integrands man- 
ifest all symmetries of the tree amplitudes. Certainly, it will be interesting to further 
study these integrands. 

Note: While this paper was being completed, works [31] appeared in which ie 
prescriptions are discussed in the context of cutting open higher-loop graphs. In the 
present work we avoided such questions by using retarded boundary conditions, which 
make the ie's clear and unambiguous, using analytic continuation in the end to convert 
to time-ordered conditions. 
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A Schwinger-Keldysh Feynman rules 

The Schwinger-Keldysh Feynman rules are summarized in Fig. 8. The propagator is 
the 2x2 matrix 

G rr G ra \ _ ( \(G> + G<) G R \ _ ( ttS(p 2 + m 2 ^ -~=± 



G ar G"J \G A j-^ 

We find convenient to draw retarded propagators, which form the two off-diagonal 
elements, with an arrow on them representing the time flow. We draw a cut on the 
(~ h) symmetrical propagator G rr , which is the only diagonal element. Interaction 
vertices with one outgoing arrow (such as rra or rrra) are the usual ones, and vertices 
with tree outgoing arrows (such as aaa or raaa) are down by Ti 2 /A. 

B The scalar triangle in D dimensions 

In this appendix we consider the triangle scalar integral in D dimensions with 3 generic 
external momenta. We want to check that the single-cut formalism reproduces it 
correctly. It is expressed in the form 



n5(q 



+ 77 . ow „ ~^7+ 



(2q-p 1 +p 2 1 )(-2q-p 2 +p 2 2 ) (2q-p 2 +p 2 2 )(-2q-p 3 +p 2 ) (2q-p 3 +p 2 3 )(-2q- Pl +p 2 

(B.58) 
We will do these integrals using Feynman parameters. The basic integral is 

7v5(q 2 ) c r 1 T(n-2+2e) 

X — —rz 7—7 — — - X 



[2q-p+k 2 ] n 2e ( p 2)i-,( A; 2)n-2+2e r(2e)r(n) ' 
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Figure 8: Schwinger-Keldysh Feynman rules, (a) retarded and cut propagators, (b) 
Ordinary vertices, that have one a index, (c) h -suppressed vertices, that have three a 
indices. 

where cr = (i^ ^va-le) • ^ ne branch of this function will be discussed momentarily. 
Proceeding formally for the moment, the first term in the bracket becomes 



/ 



«W) _ rr »w 



(2q-pi+p 2 1 ){-2q-p 2 +p 2 ) J q J [2q-{yp 1 -{l-y)p 2 ) + yp\ + {l-y)p 2 2 } 2 

cr 1 f 1 dx 



2e {pIpIY J {xp 2 + {l-x)p 2 2 - x(l-x)pj) 



2\1- 



where on the second line we have changed variables to y = {l—x)p\/{(l—x)p\ + xp 2 ). 
In Euclidean kinematics (pf > for i = 1,2,3), the denominator is real and positive 
everywhere along the naive contour which has x real. The x integral along that contour 
can be evaluated explicitly in terms of Gauss' hypergeometric function: 



("r 



x 



2e 2 (pjplpiyA^ 

l+^i Y j? , ! , l+^i x f l-iS 2 y „ ( I-182 

2 FAe, 2e, 1+e, -- — — ) - — — — 2 FAe, 2e, 1+e, 



(B.59) 



1—iSiJ 1—iSi \l+iS 2 J l+iS 2 

Here A is the is the so-called Gram determinant 

A = 2p^2 + 2p2 p 2 +2p 2 p 2 _ p A _ p A _ p A = ^2^ _ 4^.^)2 

and Si = P2 nr 3 , with similar definitions for cyclic permutations. Using hypergeometric 
identities it is possible to show that the bracket is symmetrical under interchange of pi 
and p 2 . 

The proper x contour must now be discussed. We consider, for simplicity, only the 
Euclidean regime where all external momenta are space-like p 2 > 0. The integrand has 
two branch cuts, which begin at the two complex zeroes of the denominator and reach 
out to infinity. In the second and third term in (B.58) the poles of the two factors are 
on the same side of the real axis and so the ordinary real x contour is correct. In the 
first term, the boundary condition in the two factors are not consistent which means 
that the real contour makes no sense. Rather, the contour must go from to infinity 
and then from infinity to 1, so as to ensure that, as a function of x, the pole in the 
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loop integrand never cross the real frequency axis. This contour differs from the real 
one by a discontinuity across one of the two branch cut (either one produces the same 
discontinuity). This produces the additional contribution 

dx c r i 2kT(1 - 2e) 



cr 



Summing over all contributions, we therefore find 



r 2 (i 



cr 



2e 2 (p 2 1 p 2 2 plYA 1 2- 



-[f(Si) + f(S 2 ) + f(Ss) + c] 



where 



f{Si) = - 

i 

and 



1+iSi 

1-iSi 



2 iq(e,2e,H-e, 



1+iSi 



1-iSi 



1—iSi 



1+iSi 



2 F 1 (e,2e,l+e, 



l—i5i, 
1+iSi' 



-2vre 



r(l - 2e) 



P(l-e)- 

The last term is the branch cut contribution. This agrees precisely with the form given 
in [32] , Appendix V, confirming the validity of the single-cut integrals. 

More precisely, this is the result in Euclidean kinematics. It can be converted to 
Minkowski kinematics, with either retarded or time-ordered boundary conditions, by 
means of the analytic continuations described in section 3. 

For convenience we used the same notation as [32] (up to a factor p\p\p\ in our 
definition of A). This has a finite limit as p\ and p\ go to zero, which agrees with 
(3.14), and also a finite limit as e — > which can be expressed in terms of dilogarithms 
[32]. 



C Some terms in the 2-loop vacuum polarization numerator 

We list the contributions to H^{p) with two quark lines cut. These are the forward 
amplitudes of two quarks against two currents as depicted in Fig. 7 (c). We denote q\ 
and q 2 the momenta of the two quarks. 



- / n5'(qf)n5(q 2 )Tr 



QiT 



-i u qii a Q2i a - 



qi+p' " '" ~(<?2-gi) 2 
+(G»,/i)^ (-p,u))] 



Tr5(ql)n5(q 2 2 )Tr 



1 



1 



QiT 



-7 



+qiY 



qi+p q 2 + p 

i a ! 

-7 92 7«" 



7 ?27a 



-7" 



(<?l-92) 2 

1 



qi+p' ''qi+p (qi+p-q 2 ) 2 



+ / n8(qi)*6(ti)Tr 



qiY 



1 -l a q 2 Y j 



1 



+9i7 



qi+p 
l 

t 

q 2 -p 



-la 



q 2 -p (qi+P~q 2 ) 



-'fOila 



-7 



qi+p (qi+p-q2) 2 _ 



(cm) 
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This expression is to symmetrized, with unit weight, under the obvious symmetry 
9i "H- q 2 , which we have not manifested for conciseness of notation. With this in mind, 
manifest transverseness of the integrand with respect to p^ is readily verified. 

The sum of the first and fourth line, which share the Dirac structure Tr [fcy^p-f^fp 
'-/"p-fa] upon contracting with if" , is reproduced in (5.38). The remaining contributions 
have been included in the final answer (5.40). 



D A (simple) 2-loop scalar integral 

In this Appendix we consider the 2-loop scalar integral 



91-92 



where 



Q1Q2 (<li+<l2) 2 (qi+py 



d] 



di 



r(2e)r(l - e)r(l - 2e) 

(47r)Tr(l + e)r(2-3e) 

rfe)r(i-e) 2 



(P 



2\-2e 



(D.6i; 



(4vr)Tr(2-2e) 

We want to check that this result is correctly reproduced by (2.4), which will provide 
a nontrivial check on the factor | in that equation. 

Applied to this integral, (2.4) produces the three terms 



-X" (91,92) 



91,92 



1 



+ 



+ 



(9i+92) 2 (gi+92+p) 2 (qi+q2) 2 (qi+py 
2 



(qi+q2+p) 2 (qi+py 



(D.62) 



where 



X(qi,q 2 ) 



*-\\t + {ql)S + {ql) + \t + {ql)S-{q\ 

+ ls'(qi)S + (q 2 2) + P'(q 2 1 )5-(q 2 2 ) 



The ie prescriptions in the denominators, not indicated, may be recovered trivially by 
considering the time flow in the corresponding diagrams. 

Our strategy is to first perform the integration over q 2 and then analyze the re- 
maining one over q\. The first term in the bracket depends only on the combination 
(91+92) an d for it we use that 



/ X( qi , q2)(27r) D 5 D ( qi +q 2 - q) = d T sin(Tre) 
"'91,92 



% 2 



3cos(7re) 



For the remaining two terms, we use that 



1X{q u q 2 ) 

(92+A) 2 



d r (A 



2\-t 



n5+(q 2 )+ir5-(q 2 ) _ sin^e) _ 

2 3cos7re L Wi; Wi;j 
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Actually, in the second term, A 2 = q\ = and so this whole term vanishes - - it 
contains a massless bubble with no scale in it. 

Using those results, the first and third terms of (D.62) sum to 



— dr 



d r 
' '</ 
sm^7re) 

3cos(7re, ., 



sin(7re)#(— q 2 ) %5(q 2 



q 2 \q 2 \ € (q+p) 2 (q+p) 2 ( 1+e ^ 

9{q 2 ) + 6(-q 2 ) cos(7re) . [u5 + (q 2 ) - 7r6~ (q 2 )} 



i- 



q 2 (q+p) 2 \q 2 \ 2e (g+p) 2 ( 1+e > 



(D.63) 



We have chosen to group the terms in a specific way for reasons that will soon become 
manifest. The easiest way to evaluate the integral on the first line is to note that that 
it is equal to J (— i) / q 2 (q+p) 2( - 1+e \ as follows from applying contour integration on the 
q° integral in the latter expression. In any case, it gives 

J = r(2e)r(l-e)r(l-2e) 
(47r)fr(l + e)r(2-3e) 

The integral on the second line vanishes; the easiest way to see this is to note that, in 
the frequency plane, it is equal to a horizontal contour integral that can be deformed 
to infinity without crossing any singularity. In any case, it vanishes. Equation (D.64) 
agrees precisely with (D.61), which confirms in a non-trivial way the factor | in (2.4). 

E Fast decay in M = 4 SYM 

We show that forward amplitudes in M = 4 decay like 1/q 3 at large q. The argument is 
a simple extension of the one given in [33] , to establish the absence of one- loop bubbles 
in M = 4 and supergravity. 

We begin with the representation (6.47) of the forward amplitude in terms of uni- 
tarity cuts integrands, 

A n+2 (wq) = JZ / rf V^V A L ({wq + z(I)n, T)i}, {p/ - (wq + z(I)n),r) 2 }, T)i L ) 
i J 

x /_ , ...-\-> A r({~( W( 1 + 2 ( J )^)> -Vi}, {-Pi + ( W Q + zi 1 )^), -V2}, Vi R ), 



where z(I) 



(pj+wq) 2 

2wpi-q+p% 



2pi'H 

We study each term separately. The main point is that both Al and Ar are 
themselves BCFW-deformed amplitudes: 

wq+zlljfi = w(q fi) - /i. 

Prix 2prfJ. 

The can be interpreted as a BCFW deformation, the first parenthesis giving the de- 
formation vector. Hence this places us precisely in the situation considered in [33], 
and it remains to estimate the fermion integral over 771 and r]2, or, equivalently, to sum 
over intermediate states. This was done in [33], and we briefly review their argument. 
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Using a ^-independent supersymmetry transformation, r\\ and r^ 2 can be both set to 
zero in Al and Ar at the cost of absorbing these variables into a w-independent but rj- 
dependent redefinition of the background. Al times Ar summed over a supermultiplet 
of states then become the amplitudes of four +- helicity gluons against a combination of 
supersymmetry-transformed backgrounds. Each such amplitude decays like 1/w, and, 
including the factor 1/w from the propagator, one finds A(wq) ~ 1/w 3 as claimed. 
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